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We describe a fault-tolerant one-way quantum computer on cluster states in three 
dimensions. The presented scheme uses methods of topological error correction resulting 
from a link between cluster states and surface codes. The error threshold is 1.4% for 
local depolarizing error and 0.11% for each source in an error model with preparation-, 
gate-, storage- and measurement errors. 

1 Introduction 

A quantum computer as a physical device has to cope with imperfections of its hardware. Fortu- 
nately, it turns out that arbitrary large quantum computations can be performed with arbitrary 
accuracy, provided the error level of the elementary components of the quantum computer is be- 
low a certain threshold. This is the content of the threshold theorem for quantum computation 
El E] • The threshold theorem also provides lower bounds to the error threshold which are in 
the range between 10^^'^ and 10~^, depending on the error model. It thus appears that there is a 
gap between the required and the currently available accuracy of quantum operations, and it invites 
narrowing from both the experimental and the theoretical side. In this context, significant progress 
has been made in [2] where a threshold estimate in the percent range has been demonstrated. For 
experimentally viable quantum computation there is a further desideratum besides a high error 
threshold. With the exception of certain schemes for topological quantum computation jElElEI) the 
price for fault-tolerance is an overhead in quantum resources. This overhead should be moderate. 

Here we describe a fault-tolerant version of the one-way quantum computer {QCq)- The QCq is 
a scheme for universal quantum computation by one-qubit measurements on cluster states ^9,. 
Cluster states JOl consist of qubits arranged on a two- or three-dimensional lattice and may be 
created by a nearest- neighbor Ising interaction. Thus, for the QCq, only nearest-neighboring qubits 
need interact and furthermore, only once at the beginning of the computation. For this scenario in 
three dimensions we present methods of error correction and a threshold value. 

The existence of an error threshold for the QCq has previously been established |1H I12 | [T3] and 
threshold estimates have been obtained |12| ITSj. by mapping to the circuit model. Here we take 
a different path. We make use of topological error correction capabilities that the cluster states 
naturally provide T¥ and which can be linked to surface codes [HIE]- The main design tool upon 
which we base our construction are engineered lattice defects which are topologically entangled. 

The picture is the following: quantum computation is performed on a three-dimensional cluster 
state via a temporal sequence of one-qubit measurements. The cluster lattice is subdivided into three 
regions, V , D and S. The set S comprises the 'singular' qubits which are measured in an adaptive 
basis. The quantum computation happens essentially there. The sets V and D are to distribute 
the correct quantum correlations among the qubits of S. V stands for 'vacuum', the quantum 
correlations can propagate and spread freely in V . D stands for 'defect'. Quantum correlations 
cannot penetrate the defect regions. They either end in them or wrap around them. In both cases. 




Figure 1: Topological error correction for the QCc- The three-dimensional cluster state, shown as a 
block, is structured by lattice defects which exist in two kinds: primal and dual. The defects wind 
around another. Some of them hold singular qubits (red) which realize the non-Clifford part of the 
quantum computation. 



the defects guide the quantum correlations. V and D are distinguished by the bases in which the 
respective cluster qubits are measured. The region V fills most of the cluster. Embedded in V are 
the defects (D) most of which take the shape of loops. These loops are topologically entangled with 
another. Further, there are defects in the shape of ear clips which each hold an S-qubit in their 
opening. Such a defect and the belonging S'-qubit form again a loop. These are the only locations 
where the S'-qubits occur. See Fig. ^ 

There are two codes upon which we base our construction, the planar code 15_^ and the concate- 
nated [15,3,1] quantum Reed-Muller code ^lEI- Both codes individually have their strengths 
and weaknesses, but they can be advantageously combined. The planar code has a relatively high 
error threshold of about 11% pL8| . and the symmetries of its stabilizer fit well with cluster states. 
Moreover, fault-tolerant data storage with the planar code and the creation of long range entangle- 
ment among planar code qubits are easily accomplished in three-dimensional cluster states, via a 
bcc-symmetric pattern of one-qubit measurements. Therefore, it is suggestive to base fault-tolerant 
cluster state computation on this code. However, the planar code is not well suited to non-Clifford 
operations which are essential for universal quantum computation. 

Now, in cluster- and graph state quantum computation, the non-Clifford part of the circuit 
is implemented by destructive measurements of the observables {X it Y)l\pl (The Clifford part 
is implemented by Y- and Z-measurements.). The Reed-Muller quantum code is well suited 
for fault-tolerant quantum computation via local measurements because the measurements of the 
encoded observables (Xiby)/\/2 (and of X, y, Z besides) are accomplished fault-tolerantly by the 
respective measurements on the bare level — and bare level measurements is what we are allowed 
to do in the QCq- If we could assume we had given a graph state state as algorithmic resource 
where each cluster qubit was encoded with the concatenated Reed-Muller code and where noise 
acted locally on the bare level, then fault-tolerant quantum computation were trivial to achieve. 
However, is not obvious how to create such an encoded graph state affected by local noise only. But 
this is just what the topological error correction in cluster states can do. 
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This paper is organized as follows. In Section [21 we introduce the ingredients required for the 
error correction mechanisms we use, namely cluster states, the planar code and the 15-qubit Reed- 
Muller quantum code. In Section |21 the measurement pattern used for fault-tolerant cluster state 
quantum computation is described, and in Sections IHl - El it is explained. Specifically, in Section jSl 
we describe the physical objects relevant for the discussed scheme — defects, cluster state quantum 
correlations, errors and syndrome bits — in the language of homology. In Section (HI we introduce the 
techniques for structuring quantum correlations via topological entanglement of lattice defects. Our 
error models are stated in Section 17.11 and the fault-tolerance threshold is derived in Section 17.21 
The overhead is estimated in Section El We discuss our results in Section j^l 

2 Cluster states and quantum codes 

This section is a brief review of the ingredients for the described fault-tolerant QCq- 

Cluster states. A cluster state is a stabilizer state of qubits, where each qubit occupies a site 
on a d-dimensional lattice C. Each site a G C has a neighborhood N{a) which consists of the 
lattice sites with the closest spatial distance to a. Then, the cluster state |0)c is — up to a global 
phase — uniquely defined via the generators Ka of its stabilizer 



i.e., |(/>)c = Ka\4>)c- Here, Xa and are a shorthand for the Pauli operators Gx and that 
we use throughout the paper. We refer to the generators Ka of the cluster state stabilizer as the 
elementary cluster state quantum correlations. 

In this paper, we will use as the lattice underlying the cluster state a bcc-symmetric lattice in 
tree dimensions. That is, the location of cluster qubits is given by lattice vectors 



We sub-divide the set of qubits into two subsets, the even and the odd qubits. For even (odd) qubits 
the sum of the coordinates of their respective lattice site is even (odd). 

Note that instead of with a bcc-symmetric lattice we could have equivalently started with a 
cluster state on an sc-symmetric lattice, because a cluster state on the latter is mapped to a cluster 
state on the former by Z- measurements on the qubits (e, e,e) and (o, o, o); see P|. 

The [15, 3, 1] quantum Reed-MuUer code. By this we denote a 15 qubit CSS-code based on 
the (classical) punctured Reed-Muller code 7^(1,4)* 19 . Its stabilizer generator matrix has the 
form 



Ka.= Xa (g) Zb, VaGC 



(1) 



{(o[dd],e[ven],e), (e,o,e), (e,e,o)}, 
{(e,o,o),(o,e,o),(o,o, e)}. 



odd qubits, 
even qubits. 



(2) 




(3) 



where 



Gx 



/loioioioioioioiX 

011001100110011 
000111100001111 
VOOOOOOOllllllll/ 



(4) 



and Gz is given by Gx 



Gz© (1,1, ...,1,1). 
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Figure 2: Elementary cell of the bcc lattice. Qubits live on the faces and edges of the elementary 
cell. Syndrome bits are located in the cube and on the sites. Each elementary cell has a volume of 
2 X 2 X 2 in cluster units. Black: elementary cell of the bcc lattice, blue: edges of the cluster graph. 



This code has the fairly rare property that the encoded non-Clifford gate f7z(7r/4) = exp(— ivr/SZ) 
is local ^lEl, i.e., 

15 

exp(-i^z) -(g)exp(i^Zi). (5) 

i=l 

This property has been used in magic state distillation In the computational scheme described 
here we use it to fault-tolerantly measure the encoded observables via local measurements of 
observables . 

Surface codes. For the surface codes physical qubits live on the edges of a two-dimensional 

lattice. The support of a physical error must stretch across a constant fraction (typically 1/2) of 
the lattice to cause a logical error. The protection against errors is topological. 

The stabilizer generators of the code are associated with the faces / and the vertices v of the 
lattice, 

Sx{v)= (g) Xe, Sz{f)= (g) Ze. (6) 

e\v£{de} ee{df} 

Therein, d is the boundary operator. The number of qubits that can be stored depends on the 
boundary conditions of the code lattice. The code resulting from periodic boundary conditions, the 
'toric code' jHj, can store two qubits. 

As an example we would briefly like to discuss the planar code 45 which encodes one qubit; see 
Fig. Ek- This example exhibits many features of our subsequent constructions one dimension higher 
up: Errors are identified with 1-chains and show a syndrome only at their end points. Homologically 
equivalent chains correspond to physically equivalent errors. Error chains can end in the system 
boundary without leaving a syndrome. 

Specifically, Pauli operators Zi live on the edges of the primal (=shown) lattice, and Pauli 
operators Xj live on edges of the dual lattice. The encoded Z-operator is a tensor product of 
individual Zi operators corresponding to a primal 1-chain stretching from left to right across the 
code lattice. The encoded X-operator corresponds to a 1-chain of the dual lattice that stretches 
from top to bottom. 

The code stabilizer is modified at the system boundary. For example, a face to the left or right 
of the lattice has only three elementary 1-chains in its boundary, instead of four. Such boundary is 
called a 'rough edge'. Where no modification of the faces occurs the system boundary is a 'smooth 
edge'. 'Smooth' on the primal lattice is 'rough' on the dual, and vice versa. Error chains can end 
in a rough edge of their respective lattice without leaving a syndrome, but not in a smooth edge. 
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Figure 3: Codes, a) Planar code. The encoded Pauli observables X and Z, two errors and the 
different boundary types 'smooth' and 'rough' are shown. The errors Ex, Ez end in their respective 
rough boundary. Ex corresponds to a trivial 1-cycle on the dual lattice and has no effect on the 
encoded state. Ez corresponds to an non-trivial 1-cycle on the primal lattice. It leaves no syndrome 
but causes a logical Z-error. b) Bi-colorable graph state locally unitary equivalent to the |-|-)-state 
of the 15-qubit quantum Reed-Muller code. 



The surface codes will occur rather implicitly in our constructions. The reason is that here we 
do not use such codes to encode logical qubits. Instead, we use them to appropriately "wire" a 
subset of the cluster qubits, the S'-qubits. The link between surface codes and cluster states has 
been established in Jl], for the purpose of creating long range entanglement in noisy 3D cluster 
states via local measurements. It has been found that the error correction implemented by the local 
measurements is described — like fault-tolerant data storage with the toric code — by the so-called 
Random plaquette 7j2-90.uge model in three dimensions (RPGM) 18 . The three-dimensional cluster 
state is like a surface code, one dimension higher up. The third dimension, which is temporal in data 
storage with the toric code, is spatial for the cluster state. The extra spatial dimension can be used 
to fault-tolerantly mediate interaction among qubits. The creation of an encoded Bell state over 
large distances ^1] is the simplest example. The long-range quantum correlations are engineered 
by the suitable choice of boundary conditions. 

Why are the above two codes chosen? For the fault-tolerant scheme of quantum computation 
described in this paper we require a quantum code with the following three properties: 1) The code 
is of CSS-type, 2) The code satisfies Eq. Q, and 3) The code fits with cluster states. For some 
arrangement of qubits on a translation-invariant two-dimensional lattice, the code has a translation- 
invariant set of stabilizer generators and these generators each have a small support on the lattice. 

The Reed-Muller code has properties 1 and 2 but not 3. The surface codes have properties 1 and 
3 but not 2. Thus, neither of the codes alone suffices. But their combination does, as is described 
in the subsequent sections. 

3 The measurement pattern 

As pointed out in the introduction, we subdivide the cluster C into the three disjoint subsets V, 
D and S. S is the the set of qubits where the non-Clifford part of the quantum computation is 
performed, and V and D are to connect the qubits of S in the proper way. We have not yet explained 
what the defects are, and we will do so only in the next section. For the moment it suffices to note 
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that the defects are located on the subset D of the cluster, and that D is the union of the two 
disjoint subsets Di and D2. The measurement pattern on D, S and V is given by 

Defect qubits a & D : Measurement of | , if | ^ ^ '(^^ , 

\ Za ' \ a £ Di ' 

Singular qubits a G S : Measurement of — " ^ ^'^^ 

Vacuum qubits a £ V : Measurement of Xa- 

Now we have to explain why we choose this measurement pattern, which is best done using the 
language of homology. 



4 Involving the Reed-Muller quantum code 

In this section we explain the role of the Reed-Muller code for the described computational scheme. 
Consider a cluster state \4>)c2 on a two-dimensional cluster €2- It is a resource for universal quantum 
computation by measurements of the local observables X, Z and , [S]. Denote by Q C C2 the 

set of qubits which are measured in the eigenbasis of . These measurements implement the 

v2 

non-trivial part of a quantum circuit. The measurements of X and Z on the qubits C2\Q implement 
the Clifford part. They are performed simultaneously in the first round of measurements. |^aigo)Q 
is the state of the unmeasured qubits after the first measurement round. It is an algorithm-specific 
stabilizer state, hence the subscript "algo". Since it is a stabilizer state, it is easy to create and 
one may start with this state as an algorithm-specific resource instead of the universal cluster state. 
Quantum computation with this state proceeds by measuring local observables . 

Now suppose an encoded version of this state, |^aigo)s! was given. The state were not perfect 
but only affected by local noise on the bare level. Of course, a question that arises immediately is 
how such a state is obtained. The main part of work in this paper goes into answering this ques- 
tion, see subsequent sections. Now, with |^'aigo)5 given, one could perform fault-tolerant quantum 
computation by fault-tolerant measurement of the encoded observables . This is not what 
we have in mind, because we are seeking a scheme of fault-tolerant quantum computation by local 
measurements. Here lies the reason for involving the (concatenated) Reed-Muller code: For this 
code, the fault-tolerant measurement of the observables proceeds by local measurements of 

observables . The reason for this is property ((Sj). If J is a set such that gx{J) '■= ^jej-^j is 
in the RM code stabilizer, then also 

X ■ ^Y- 

9±{J) ■■= 6^ ^ r- ^ ^ code stabilizer. (8) 



The relevant encoded observables are given by 

X ± F ^ X. =F Yj 



V2 ^ \/2 



(9) 



The "=f" is for the total number of concatenation levels being odd. If the number is even, replace 
"=F" by "lb". Therefore, if all the bare qubits belonging to an encoded qubit are individually 
measured in the eigenbasis of ^^^^ , then the eigenvalue found in a measurement of the encoded 

observable and the eigenvalues of the stabilizer generators g+{J) can be deduced from the 

individual measurement outcomes. This is all what is needed for fault-tolerant measurement of the 
encoded observable ■ 
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The Z-part of the code stabihzer is lost in the local measurement, but it is not needed for the 
fault-tolerant measurement of . This can be seen as follows. For simplicity assume a local 

depolarizing error p/3 {[Xj] + [Yj] + [Zj]) = ^ {[Xj + Yj] + [Xj - Yj] + V2[Zj]) for all qubits j. 
The brackets "[•]" indicate a super-operator. W.l.o.g. assume that the local measurements are in 
the eigenbasis of Then, the error ^^^^ is absorbed in the measurement and has no effect. 

The second error ^^^^ = i ^^^^ — ^j: so all remaining errors are equivalent to Z-errors. Such 
errors are identified by the stabilizer elements {5+ (J)}. 

To summarize, the Reed-Muller code is involved to perform the fault-tolerant measurement of 
the encoded observables locally on the bare level. The eigenvalues corresponding to the encoded ob- 
servable and to the relevant stabilizer generators are simultaneously inferred from the measurement 
outcomes. Error correction proceeds by classical post-processing of these quantities. 

5 Involving a topological quantum code - Homology 

The remaining question is how we actually create the state l^aigo)^ with only local error from a 
three-dimensional cluster state. To accomplish this task we involve topological error correction. 

5.1 Errors and correlations as chains 

The physical objects of our discussion — cluster state correlations and error operators — may be 
identified with faces and edges of an underlying lattice. Compositions of such faces or edges are 
called 2-chains and 1-chains, respectively. For the chains homology provides an equivalence relation; 
namely, two chains are homologically equivalent if they differ by the boundary of a third chain one 
dimension higher up j2Ul 1211 1^ . Homology plays a role in our constructions because homological 
equivalence of the underlying chains implies physical equivalence of the associated physical objects. 
First we introduce the two simple cubic sub-lattices C and C whose vertices are at locations 

r ■ 11' (10) 

C : {(0,0,0)}. 

One lattice can be obtained from the other via translation by a vector (ibl,ibl,ibl). £ and C are 
dual to another in the sense that the faces of C are the edges of £, the cubes of C are the vertices 
of C, and vice versa. We denote by * the duality transformation that maps primal edges into the 
corresponding dual faces (*e = /), and so forth. 

We denote by B{Co) := {vk} the set of vertices in £, by I3{Ci) = {ei} the set of edges in £, by 
S{C2) = {fm} the set of faces in C, and by 8(03) = the set of elementary cells [or cubes] of 
jC. We may now define chains in C [20]. l3{Co) forms a basis for the set Cq of so called 0-chains cq, 
B{Ci) forms a basis for the set Ci of 1-chains ci, and so forth. Specifically, the chains are given by 

Co = T^kZkVk, Ci = YliZieu C2 = J2rn^mfm, C3 = En ^n9n- (H) 

where Zk, zi, Zm, Zn G ^2- The sets Co, Ci, C2 and C2 are, in fact, abelian groups under component- 
wise addition, e.g. ci + c'l = Yli^i^i + Yli^'i^i — Ylii^i + Fo'^ each i = 1..3, there exists a 
homomorphism di mapping Ci to Cj_i, with the composition o 5^ = 0. Then, 

C = {C3,C2,Ci,Co] (12) 

is called a chain complex, and d is called boundary operator. It maps an i-chain q to its bound- 
ary, which is an i — 1-chain. In the same way, C can be defined as a dual chain complex, 
C = {C3, C2, Ci, Co}, with chains C3, C2, ci and cq. 
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Now, considering a space C, two chains c^, G Cn{C) are homologically equivalent if c'„ = 
c„ + 5c„,+i for some c„+i G C„,+i(C) 20 . Of interest for topological error-correction is the notion 
of relative homology. Consider a pair of spaces {C,D) with D C C. Then, two chains 
called equivalent w.r.t relative homology, =1- c„, if = Cn + dcn+i + 7n for some c„+i G C„+i(C), 
7„ G Cn{D); see j^J- Relative cycles may end in D. 

Below we describe how the cluster state quantum correlations may be identified with the 2- 
chains, the errors with the 1-chains and the syndrome with the 0-chains of C and C. All these 
objects appear in two kinds, 'primal' and 'dual', depending on whether they are defined with 
respect to £ or £. 

Cluster state correlations. We define primal such correlations, K{c2), and dual ones, K{c2), 
which can be identified with 2-chains in C and C, respectively, by 

Kic2):= n K{C2):=JI Kj. (13) 

/G{c2} 7e{c2} 

Therein, e.g. the set {02} is defined via a mapping 02 = Y^m^rnfm — > {C2} = {fm\zm = !}• 
Further, we introduce the notion 0(c) := (S>ae{c} ^a, for all c G C and O G {X, Z}. It is now easily 
verified that 

K{C2) = X{c2)Z{dc2), K{C2) = X{c2)Z{dc2). (14) 

Errors. We will mainly discuss (correlated) probabilistic noise. Then, it is sufficient to restrict 
the attention to Pauli phase flips Z, because Xa = 'S)beN{a) etc. We combine Z-errors on odd 
(even) qubits to primal (dual) error chains E{ci) {E{ci)), 

E{ci):=Z{ci), E{ci):=Z{ci). (15) 

Syndrome. The first type of correlations we discuss are those for error correction in V . They are 
characterized by the property that the corresponding 2-chains have no boundary. I.e., we consider 
-fC(c2), Kic2) with {c2},{c2} G V and dc2 = dc2 = 0. With H13() . these correlations take the form 
K{c2) = ^(c2), K(^2) = X{c2)- They are measured by the X-measurements in V, see Q, and are 
used to identify errors occurring on the qubits {02}, {02}- 

The group of 2-chains in the kernel of d we denote by Z2{C). Since 95 = 0, a subgroup of those, 
-62(^2), is formed by the 2-chains which are themselves a boundary of a 3-chain. Denote by q (q) a 
3-chain from the basis ^(Cs) {B{Cs)). It represents an individual cell [or cube] of the primal lattice 
£ (dual lattice £). The associated quantum correlations are 

K,:=K{dq)=X{dq), := K (Oq) = X (dq) . (16) 

When being measured, each of these correlations yields a syndrome bit Sy{q), Sy{q) which, we say, 
is located at q or q, respectively. Because the lattices C and £ are dual to another, we may identify 
the cell g in £ with a vertex v in £, and vice versa. In this way, the syndrome bits become located 
at vertices of the lattices £ and £. The syndrome resulting from the quantum correlations (|16|1 are 
those which enable topological error correction |18j . 

Syndrome and errors. Let E{ci) denote a primal error chain, q G B{C3) a cell in the dual 
lattice £ and v G 13(Co) = *q- K{q) detects the error E{ci) if |{ci} n = odd. Equivalently, 

Kiji) detects E{ci) if v G {dci}. Thus, error chains show a syndrome only at their ends. 
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Correlations and errors. Primal cluster state correlations are affected by dual error chains and 
dual correlations are affected by primal error chains. Primal correlations are not affected by primal 
error chains, and dual correlations are not affected by dual error chains. 

To see this, note that a primal correlation K(c2) consists of Pauli operators X on even qubits 
and Pauli operators Z on odd qubits. A dual error chain E{ci) consists of operators Z on even 
qubits. Then, E{ci)K{c2) = (-l)l{Si}u{c2}| K{c2)E{ci). If |{ci}U{c2}| is odd, the correlation K{c2) 
is conjugated to —K{c2) by the error. If it is even, then the correlation remains unchanged. This 
situation has a geometric interpretation. |{ci} U {c2}| is the number of intersection points between 
the primal 2-chain C2 and the dual 1-chain ci. If the number of intersections is odd (even) then the 
correlation is (is not) affected by the error. For dual correlations and primal errors the situation 
is the same. Further, a primal error chain consists of Pauli operators Z on odd qubits. Thus, 
[K{c2),E{ci)] = always. Similarly, [K{c2), E{ci)] = 0, Vci,C2. 

Defects. The purpose of defects is to structure the space underlying the pair of lattices C,C 
Practically, a defect can be thought of as a set of qubits that are removed from the initial cluster 
C before the remaining qubits are entangled. For the chain complexes C, C, a defect is a set d of 
missing edges. What defines a defect as an entity is that the belonging edges are connected. As all 
the other objects, defects are either primal or dual, 

d C B{Ci), d c B(Ci). (17) 

The sets d, d of defect qubits are not arbitrary. Seen from afar they take the shape of doughnuts. 
These doughnut-shaped defects will be topologically entangled with another, and the way they 
are entangled encodes the quantum algorithm to be performed. From the viewpoint of quantum 
logic, what matters about the doughnuts is that they are loops. Their 'thickness' is required for 
fault-tolerance. 

We now briefly explain how the above definition of a defect as a set of missing cluster qubits fits 
with the measurement pattern 0. Formally, each defect d will be assigned a set D{d) of locations 
on the cluster. This set is subdivided into a set of edge- and a set of face qubits, Di{d) and D2{d). 
Here, the notions of 'edge' and 'face' are in reference to the lattice the defect belongs to. If the 
defect is primal (dual) then the edges and faces are taken with respect to the primal (dual) lattice. 
For primal defects, the sets Di[d) and D2{d) are defined as 

D^{d) := d, D2{d) = {/ G B{C2)\ {df} nd = {df}} . (18) 

For dual defects, replace / by / and C2 by C2 in the above definition. The whole defect region D 
splits into an edge part Di and a face part D2 , D = Di U D2 where 

D^ = [jD^{d), D2 = ]jD2{d). (19) 
d d 

Now the measurement pattern becomes understandable: the edge qubits in the defects are 
measured in the Z-basis which effectively removes them from the cluster ^Uj. In this way, the 
quantum state on the exterior of the defect becomes disentangled from the state with support on 
interior of the defect. Thereby, a defect is created in the cluster lattice. Note that the qubits on 
faces whose entire boundary is in Di{d) become disentangled individually. If no errors were present 
we could leave these qubits alone. However, their measurement in the X-basis provides additional 
syndrome and so it is advantageous to measure them. 

Correlations and defects. In the proximity of a primal defect, edges in the boundary of a primal 
2-chain are removed, see (|17|). Therefore, primal correlations can end in primal defects. Dual defects 
do not remove primal edges, and thus primal correlations cannot end in dual defects. Analogously, 
dual correlations can end in dual defects, but not in primal defects. 
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to this onej. 


dual corr. 


primal defect 


This object do 

dual defect 


es ... 

primal err. cy. 


dual err. cy. 


correlation 


nothing 


bound 


repel 


nothing 


affect 


dual correlation 




repel 


bound 


affect 


nothing 


primal defect 






encircle 


pairwise end 


encircle 


dual defect 








encircle 


pairwise end 


primal err. cy. 










nothing 



Table 1: This table shows who does what to whom. 'A bounds i?' is synonymous with 'i? ends in 
A\ The displayed objects do not interact with themselves. 

Syndrome and defects. In the presence of a primal defect d the correlations {Kg = K{dq) \ q G 
B{c3) A dq n d ^ (/>} do not commute with the measurements 0, such that the syndrome at the 
locations 

Do = {v e de\ e G 4 (20) 

is lost. Note, however, that for each defect d there will be one syndrome bit associated with 
the defect as a whole. There exists a 2-cycle C2{d), {c2{d)} C V, that wraps around d, and 
K{c2{d)) = X{c2{d)). When the qubits in V are measured in the X-basis, this correlation yields an 
additional syndrome bit. Dual defects act analogously on the dual lattice. 

Errors and defects. Because the local syndrome is lost at the surface of a defect, primal error 
chains can potentially end in primal defects. However, there is a dual correlation K(c2{d)) = 
X(c2id)) wrapping around a primal defect, and this correlation detects a primal error chain E{ci) 
if the number of intersection points between C2{d) and ci is odd. Thus, primal error chains can 
pairwise end in primal defects. 

Primal error chains cannot end in dual defects, because dual defects do not remove primal 
syndrome. Similarly, dual error chains can pairwise end in dual defects, and they cannot end in 
primal defects. 

The relations among cluster state correlations, errors and defects are summarized in Tab. Q 
5.2 Homological and physical equivalence 

We have so far identified physical objects — correlations and errors — with chains of a chain complex. 
In this section we point out that it is the homology class of the chain rather than the chain itself 
which characterizes the respective physical object. The equivalence of two chains under relative 
homology implies the physical equivalence of the corresponding physical operators. 

1. Cluster state correlations. We regard two cluster state correlations K{c2), K^c^) as physically 
equivalent if they yield the same stabilizer element for the state |^'aigo)s after the measurement of 
the qubits in V and D. This requires two things. First, K{c2) and -/^(cg) need to be simultaneously 
measurable. With Oa the locally measured observables ((7| we require [K{c2),Oa] = Va G F U 
D <^=^ [i^(c2), Oa\ = Va G VUD. Second, the two operators must agree on S, K{c2)\s = -^(c2)|s- 
Then, the following statement holds: If c'2 =r C2 w.r.t. {VVJ D,D) then K{c'2) = K{c2). 
Proof-. There exists C3 G C3I {dc^} dV U D2 and 72 G C2I {72} C D2 such that c'2 = C2 + 9c3 + 72. 
1. Simuhaneous measurability: K{-^2) = {^aeD2{i2) (*^fceZ)i{72) -^^')' where 1)1(72) C Di and 
-^2(72) C D2. Therefore, with ©, [K (72), Oa] = Va G F U L> (*). Similarly, K{dc^) = X{dc^) 
such that, with ©, [K{dc^), O^] = Va G F U L» (**). Since Kic'^) = K{c2)K{dc3)K{-f2), (*) and 
(**) imply simultaneous measurability of K{c2) and K{c'2) on V ^ D. 2. Same restriction to S: 
K{dc3) and i^(72) don't act on S, hence K{c2)\s = K{c2)\s- D 
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2. Errors. Two errors E{ci) and E{c'i) are physically equivalent if they cause the same damage 
to the computation. That is, they have the same logical effect and leave the same syndrome. Then, 
the following statement holds: If c[ ci w.r.t. {V U D,D) then E{c'^) ^ E{ci). 
Proof: There exist C2 € C2, 71 G Ci, with {02} C ^ U 1)2, {71} C .Di, such that c'^ = ci + dc2 + 71- 
Now, a Pauli spin flip error X is absorbed in a subsequent X-measurement and has no effect on the 
computation, = ±^I. Thus, with Q, Xa = for ah a G F U L'2. Similarly, = h for 

all 6 G Then, E{ci + 9c2 + 71) = E{ci)Z{dc2)Z{-ii) = E {ci)K (02) X (02) Z{-fi) ^ E{ci). □ 

6 Constructive techniques 

The purpose of the measurements in V and D is to create on S the Reed-Muller-encoded algorithm- 
specific resource |^aigo)5 described in Section |1J In Section we specify the location of S-qubits 
with respect to the lattice defects and then, in Section [6.21 we give a construction for a topologically 
protected circuit providing |\I'aigo)s- 

6.1 Location of the S'-qubits 

The S'-qubits have very particular locations within the cluster. Besides the defects in the shape of 
doughnuts that we have already introduced the cluster also supports defects shaped like ear clips. 
The opening of these ear clip defects is only one cluster qubit wide. If this one cluster qubit were 
a defect qubit [q G D) too, the ear clips would become doughnuts. But the particular cluster qubit 
is not in D, it is an S-qubit. The situation is displayed graphically in Fig. 

The appropriate stabilizer generators among the 5-qubits are induced from the cluster state 
correlations associated with relative 2-cycles, by measurement of the V- and D-qubits. As everything 
else in this computational scheme, the 5-qubits occur in the two kinds 'primal' and 'dual'. We call 
an S-qubit q primal, q G Sp, if it lives on the a face of the primal lattice, and we call it dual, q G Sd, 
if it lives on a face of the dual lattice. 

We now discuss how primal and dual correlations affect the S'-qubits. Consider, for example the 
correlation K{c2) corresponding to a primal relative 2-chain C2- A primal S-qubit at location q may 
lie within C2, but never in its boundary, {802} n g = 0. A dual S-qubit q' may lie in the boundary 
of a primal 2-chain C2 but never in C2 itself, {02} H = 0. We therefore conclude that 

A primal correlation K{c2) acts on a primal S-qubit by one of the two Pauli- operators 
X, I and on a dual S-qubit by one of the two operators Z , I. A dual correlation 
K(c2) acts on a primal S-qubit by one of the two Pauli- operators Z, I and on a dual 
S-qubit by one of the two operators X, I. 

For finding the extended relative 2-cycles on the primal lattice, we thus regard Sp as part of V 
and Sd as part of D. Analogously, for finding the extended relative 2-cycles on the dual lattice, we 
regard Sp as part of D and Sd as part olV. In this way, the problem of finding the extended primal 
2-cycles on a cluster with S-qubits is reduced to the same problem without S-qubits. 

6.2 Creating l^'aigo) among the S'-qubits 

The construction of a topologically protected circuit providing |^'aigo)s' proceeds in three steps. 
First we show that |^'aigo)5 is local unitarily equivalent to a bi-colorable graph state, by local 
Hadamard-transformations. Second, we show how to create a bi-colorable graph state. Third, we 
take care of the Hadamard-transformations. 

Equivalence of l^'aigo) to a bi-colorable graph state. Denote by \^)c2 the state after the first 
measurement round {X- and Z- measurements) in a QCc-computation on a cluster state \4>)c2j cf. 
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Figure 4: a) Location of a singular qubit (dual, held by a primal defect). b,c) The correlations K, 
K affect the 5-qubit by a Pauli-operator X in (b), and by Z in (c). 



Section^ Denote by Cx, Cz the subsets of C2\Q whose qubits are measured in the X- or Z-basis, 
respectively. Further, denote by Cevem Codd the sets of even and odd qubits in C2 (checkerboard 
pattern). The state \^)c2 is given by 



(22) 



and has the form \^)c2 = l^aigo)(3 ® kest)c2\Q- Further, \(l))c2 is local unitary equivalent to some 
CSS-state, |(/.)c2 = {<^^ec,^^ Hi) \CSS)c2. Then, with ^ and C'^ := {Cx n Ccvcn) U {Cz n Codd), 




C' 



z - {Cz n c, 



u {Cx n c, 



odd J 



1^; 



C2 



I±Xa 







--even 



(^•g^ dd ^0 |C'5'5')c2- Thus, also the state I^I'aigo) is l.u. equivalent to a CSS-state 



\CSS) 



C2 



1^ 



algo/ 



Ha \CSS") 



(23) 



v<?6CoddnQ 



Now, we consider the concatenated-Reed-MuUer-encoded resource |^'aigo)5i which may obtained 
from the bare state |*aigo)Q via encoding, |^'aigo)5 = (^gegEnCg^ |^aigo)Q- The encoding pro- 
cedure Enc takes every qubit g e Q to a set S{q) of qubits, UgGQ '^(^) ~ ^- property 
that Encg o Hq = {^i(zs{q) Hi) o EnCg. Therein, Enc' is an encoding procedure for the code con- 
jugated to the Reed-Muller-code, i.e., for the code with the X- and the Z-block of the stabilizer 
interchanged. This code is of CSS type, like the Reed-Muller code itself. The encoding procedure 
changes when passing through the Hadamard-gate because the encoded Hadamard-gate is not local 
for the Reed-Muller quantum code. 

At any rate, the state |^aigo)s is l.u. equivalent to a CSS-state encoded with CSS-codes, i.e., to a 
larger CSS-state, |^'aigo)5 = {'S>qeCoddnQ '^ieS{q) Hi) \CSS"')s- Every CSS-state is l.u. equivalent 
to a bi-colorable graph state ^23 , by a set of local Hadamard-transformations. Thus, we finally 
obtain 

|^algo)5 = (g) \r)s. (24) 

Therein, Sh is some subset of S and \T)s is a bi-colorable graph state with adjacency matrix of the 
corresponding graph 
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primal dual 




Figure 5: Creation of an arbitrary bi-colorable graph state. The singular qubits are displayed as 
red circles. G denotes the sub-matrix of the adjacency matrix T which encodes the edges between 
primal and dual vertices. 



Circuit for a bi-colorable graph state. The circuit layout for the topologically protected 
creation of an arbitrary bi-colorable graph state is shown in Fig. [21 The circuit consists of a set 
of horizontal primal and a set vertical dual defects. Each primal defect comes close to each dual 
defect once, and the two defects may be linked in that region. The form of this junction is decided 
by a corresponding element of the graph state adjacency matrix G. If Gij = 1 then the defects are 
linked and otherwise they are not. In addition, each of the loop defects is linked with an ear-clip 
shaped defect which holds an 5-qubit. The graph state in question is formed among these qubits 
after the remaining qubits have been measured. 

An explanation of the circuit in Fig. El for the specific example of a line graph, is given in Fig. (HI 
From this example it should be clear how the circuit works in general. For the line graph we have 
the adjacency sub-matrix 

/ 1 1 \ b 
Giine = \^ 1 1 y d ■ (26) 

ace 

This implies, for example, that in the circuit of Fig. the dual defect winding around the dual 
S'-qubit a will be linked with the primal defect winding around the primal S'-qubit b. Now we 
explain how the stabilizer element Kb = ZaX^Zc for the graph state iFune) emerges. The other 
stabilizer generators emerge in the same way. 

Consider the relative 2-cycle 02(6) and imagine it being built up step by step. We start around 
the 5-qubit b. Because C2{b) is primal and b is primal, K{c2{b)) affects qubit & by a Pauli-operator 
X, see Eq. (PT|l. Also, 02(6) is bounded by the primal defect encircling b. 

We move further to the left. At some point, 02(6) approaches a dual defect. Primal correlations 
cannot end in dual defects. Therefore, 02(6) bulges out and forms a tube wrapping around the 
dual defect, leading downwards. It ends in the primal defect holding the dual 5-qubit c. K{c2{b)) 
affects qubit c by a Pauli-operator Z. Back at the junction, 02(6) continues to expand to the left. 
It approaches a second dual defect where it forms another tube. In result, the dual 5-qubit a is 
affected by a Pauli-operator Z. Further to the left, the primal defect closes up and bounds K{c2{b)). 

K[c2(b)) takes the form Xj^Z^Zc (^^g^/^^-jy for some set V{b) C V. All qubits in V are 
measured in the X-basis ^ such that after these measurements the correlation = zizZaX^Zc 
remains. This is a stabilizer generator for the graph state in Fig. since (a, b) and (c, b) are the 
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a) b) 




Figure 6: Explanation of the circuit in Fig.[Slfor a particular example, the line graph, a) The relative 
2-cycle 02(6) which gives rise to the graph-state stabilizer Ki, = ZaXi,Zc- b) Graph corresponding 
to the created graph state. 



only edges of the line graph ending in the vertex b. For every edge in the graph there is a link 
among defect loops in the circuit. 

The proof for the programmable circuit of Fig. [S] realizing a general bi-colorable graph state is 
similar. As an outline, each stabilizer generator is associated with a doughnut-shaped defect in the 
circuit. Such a defect bounds a correlation, and this correlation affects one 5-qubit by X. Further, 
because the considered defect is linked with other defects of the opposite kind, the correlation 
surface forms tubes. These tubes affect one other S'-qubit each — the neighbors of the first — by 
Pauli-operators Z. 

Implementing the local Hadamard-transformations. The equivalence between the graph 
state ir)^ and the encoded algorithm specific resource |^aigo)5 is by Hadamard-transformations 
on some subset Sh of S'-qubits. Wherever such a Hadamard-transformation needs to be applied, 
attach an extra loop to the circuit in Fig. El 




(27) 



7 Error sources, error correction and fault-tolerance threshold 
7.1 Two error models 

Below we describe two error models, a simple one and a more realistic one. 

Error Model 1. The cluster state is created only on those cluster qubits which are needed for the 
computation, i.e., onV \J S. The defect qubits of D are left out. 

1. The noise is described by independent partially depolarizing channels acting on each cluster 
qubit. The noisy state pc is given by pc = ^^Ta{pi) |(/>)c((/'|, with 

aec 

TPip,) = {l-pi)[Ia] + I {[Xa] + [Ya] + [Z,]) . (28) 

2. The classical computation for syndrome processing is instantaneous. 
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The reason for considering this error model first is its simphcity. We would like to separate the 
intricacies inherent in the presented error-correction scheme from additional difficulties incurred by 
a realistic error model. The basic justification for such an approach is this: The error correction 
used here is topological. Therefore, a threshold should exist regardless of whether independent 
errors are strictly local or only local in the sense of having a support of bounded size. 

The most straightforward method to create a cluster state is from a product state |+)a via 
a (constant depth) sequence of A(Z)-gates ^Uj. If these gates are erroneous, then Error Model 1 
does not apply in general.^ Specifically, one may raise the following objections to Error Model 1: 

• No correlated errors are included in Error Model 1. Creating the cluster state via a sequence 
of gates will, however, lead to correlated errors in the output cluster state. 

• Storage errors accumulate in time. There is temporal order among the measurements such 
that the computation takes a certain time tcomp which cannot be bounded by a constant for 
all possible computations. As a consequence, for the qubits measured in the final round the 
local noise rate increases monotonically with tcomp and exceeds the error threshold. 

• To leave the D-qubits out is a deviation from the originally envisioned setting: the cluster 
state on U y is algorithm-specific.^ 

To account for these inadequacies, we consider a second error model. 

Error Model 2. A cluster state on a bcc- symmetric lattice is created in four steps of nearest- 
neighbor A{Z)-gates. The gate sequence is as shown in Fig. ^ Errors occur due to the erroneous 
preparation of the initial \+) -qubits, erroneous A{Z)-gates in the process of creating the cluster state, 
storage and measurement. 

1. The computation is split up into steps which performed on sub-clusters Ck- In each step, 
unmeasured qubits remaining from the previous step — the hand-over qubits — are loaded into a 
cluster state on a sub-cluster. Subsequently, all but a few cluster qubits (the new hand-over 
qubits) are measured. The steps have their temporal depth adjusted such that each qubit, after 
being locally prepared and entangled, waits at most a constant number to of time steps until 
its measurement occurs, to > 1. Error in storage is described by a partially depolarizing noise 
with error probability ps per time step. 

2. The erroneous preparation of initial \+) -qubits is modeled by the perfect procedure followed 
by local depolarizing noise (I28() . with probability pp. Measurement is described by perfect 
measurement preceded by partially depolarizing noise with error probability pM- The erroneous 
A{Z)-gates are modeled by the perfect gate followed by a 2-qubit depolarizing channel 

T^}{P2) = (1 - P2)[Iej] + ff (lie ® Xj] + ... + [Z^ Zj]) . (29) 

3. Classical syndrome processing is instantaneous. 

^There may be situations in which the error Model 1 is in fact a good approximation. For example, consider 
a scenario in which the cluster state is purified before being measured for computation. Of course, the gates in a 
purification protocol would be erroneous, too, such that the purified state is not perfect. In effect, the errors of the 
initial state were replaced by the errors of the purification protocol. There exist purification protocols |24| in which 
the gates act transversally on two copies of the cluster state (one of which is subsequently measured). As a result, 
the errors introduced by the purification are approximately local, as in Error Model 1. The purification protocol I24| 
in its current form has a problem of its own, though; due to the exponentially decreasing efficiency of post-selection, 
it is not scalable in the size of the state. But chances are that this can be repaired. 

^For the creation of the cluster state this makes little or no difference: fewer gate operations are needed than for 
the creation of \(f>)c- For parallelized procedures that make use of the translation invariance of C, it should not be too 
difficult to remove the superfiuous D-qubits from the lattice before the remaining qubits are entangled. 
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Step 1 Step 2 Step 3 Step 4 



Figure 7: Steps of A(Z)-gates for the creation of a cluster state on a bcc-symmetric lattice, for the 
sub-clusters Ck- For k odd, the sequence is 1 — > 2 — > 3 ^ 4, and for k even it is 3 ^ 4 ^ 1 — > 2. 



In the subsequent sections we compute a fault-tolerance threshold for both error models. In 
Error Model 2, storage error will cause the minimum damage for the smallest possible value of to, 
which is io = 1- In this case, each sub-cluster carries a subset of S-qubits with no mutual temporal 
dependence. All qubits in V and D, except the hand-over qubits (See Appendix^, are measured 
immediately after being entangled. The measurement of the S'-qubits has to wait one time step. 
Henceforth we set to = 1- 

7.2 Methods for error correction and the threshold value 

There are different methods of error correction associated with the different regions V, S and D of 
the cluster. In D, there are as many inequivalent errors as there are syndrome bits, such that the 
error correction is trivial. The error correction in V is based on the random plaquette Z2-gauge 
model in three dimensions ^Hl- The error correction in S is carried out using the (concatenated) 
quantum Reed-Muller code. 

Error correction in D. In the domain D{d) of the defect only Z-errors matter for face qubits 
and only X-errors matter for edge qubits. Any other errors may be absorbed into the subsequent 
measurements ((T)). The X-errors on edges may be relocated to Z errors on the neighboring face 
qubits via = K^Xe, such that we need to consider Z-errors on face qubits only. For each face in 
d we learn one syndrome bit, yielding a unique syndrome for each error configuration. 

However, if an error X^ on an edge qubit e in the surface of defect d is relocated, the equivalent 
error K^Xg may partially be outside D[d). After error correction in D{d), an individual Z-error on 
a face qubit in V is left next to D{d). This affects the error correction in V near D; see below. 

Error correction in V. First consider a scenario where the entire cluster consists of the region 
V (i.e., there are no defects and no singular qubits). Error correction on the primal lattice C and 
the dual lattice C run separately. Here we consider error correction on the primal lattice only; error 
correction on the dual lattice is analogous. 

The error chains live on the edges of the lattice C and leave a syndrome at the end points, which 
are vertices of C. This is exactly the scenario which has been considered for topological quantum 
memory in ^Hj, and subsequently the results of JS| that we need in the present context are briefly 
summarized. The connection between topological error correction and cluster states has been made 
in ^1] for the purpose of creating long-range entanglement in the presence of noise. 

Given a particular syndrome and an error chain E(ci) compatible with this syndrome, we are 
interested in the total probability P(ci) of the homology class of ci, 

^ci) = + (30) 

where p{c[) is the probability of an individual error chain E{c'i), and the sum is over all 1-cycles. 
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For error correction we infer that the physical error which occurred was from the homology class 
with the largest probability. 

If the errors on the lattice edges occur independently with a probability q then the problem of 
computing P{ci) for a given chain ci can be mapped onto a problem from statistical mechanics, 
namely the random plaquette Z2-gauge model in three dimensions JS|- The crossover from high 
fidelity error correction at small error rates to low fidelity error correction at high error rates 
corresponds to a phase transition in this model. A numerical estimate of the critical error rate is 
qc = 0.033 ±0.001 

As far as is known, the classical operational resources required to find the most likely error 
homology class consistent with a given syndrome scales exponentially in the number of error lo- 
cations. The assumption of the classical processing being instantaneous cannot be justified under 
these conditions. However, it is possible to trade threshold value for efficiency in the error correction 
procedure. A reasonable approximation to the maximum probability for a homology class of errors 
is the probability of the lowest weight admissible chain. The minimum-weight perfect matching 
algorithm [261 127j computes this chain using only polynomial operational resources. A numerical 
estimate to the threshold with this algorithm for error correction is q'^ = 0.0293 it 0.0002 128 . 

Remark: The topological error threshold is estimated in numerical simulations of finite-size 
systems. For this purpose, the probability of logical error is plotted vs. the physical error parameter 
for various system sizes. For sufficiently large lattices (such that finite-size effects are small), we 
expect these curves to follow a universal scaling ansatz near the threshold such that they share a 
common intersection point and their slopes are proportional to a common power of the lattice size. 
As the system size is increased to infinity, we then expect the curves to approach a step function 
which transitions at the threshold value of the physical error parameter. 

The above quoted threshold value is for independent errors on the edges of £. Do the models 
for the physical error sources of Section 17.11 lead to such independent errors? The answer is 'y^s' 
for error Model 1 and 'no' for error Model 2. For the latter, we need to consider a modified RPGM 
with correlated errors among next-to-nearest neighbors. Specifically, for Error Model 1 the relation 
between the local rate pi of the physical depolarizing error and the error parameter q that shows 
up in the RPGM is 

2 

q = —pi, for Error Model 1. (31) 
3 

Given a threshold of g = 2.93% |28| for error correction via the minimum weight perfect matching 
algorithm in the bulk, then the corresponding depolarizing error rate that can be tolerated is 

pi^^ = 4.4 X 10"^ {inV). (32) 

For Error Model 2, first consider the case where only the A(Z)-gates are erroneous, pp = Ps = 
Pm = 0. Then, in addition to local errors with a rate qi there exist correlated errors with error rate 
q2 for each pair of opposite edges in all faces of C. That is, with a probability q2 simultaneous errors 
are introduced on opposite edges of the faces in C. The local noise specified by qi and the two-local 
noise specified by q2 are independent processes. The relations between the error parameter p2 of 
the A(Z)-gates and the parameters qi,q2 of the RPGM with correlated errors are 

32 / 8 \ / 64 2 / ^ V\ 

l5^^ V - T5^V 1^^^^ + I' - l5^V ) ' (33) 
1 A 4 4 ^, 2x 

The correlation of errors on sites separated by a distance of two arises through error propagation in 
the creation of the cluster state. Correlations among errors on next-neighboring sites play no role 
because such errors live on different lattices {£, and C) and are corrected independently. 
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Figure 8: Threshold estimation in lattices of finite size, for periodic boundary conditions. Here, 
with A(Z')-gates as the only error source, we find a threshold of p2,c = 9.6 x 10~^. In the inset, 
best fits to the universal scaling ansatz are drawn. Error bars denote two standard deviations due 
to finite sampling size. 



The only effect of pp,ps,PM > is an enhanced local error rate qi. We give the relations to 
leading order only; they read 

32 2 
qi = + i;{pp+Ps+ Pm) , 

^ for error Model 2. (34) 

q2 = —P2, 
15 

See Fig.lHlfor a simulation of error correction under faulty A(Z)-gates as the only error source, which 
gives rise to correlated noise between neighboring edge qubits. If we define x = {p2 —P2,c 

IlV'^o then 

the universal scaling ansatz states that fidelity F should be a function dependent only on the scaling 
parameter x in the vicinity of the threshold ^Hl. We find very good agreement (with > .9991) 
for F = A + Bx + Cx^ , where we fit for constants A, B, C, P2,ci and vq. This gives very tight bounds 
on the critical probability p2,c = 9.6 x 10~^. Interestingly, we also find = 1.00 it 0.02, which 
indicates that this model belongs to the same universality class as the purely local error model of 
the 3D-RPGM [HI- 

In Fig. inithe threshold trade-off curve between pp +ps +Pm and p2 is displayed. Numerically, 
we obtain for the thresholds in the bulk 

Pp = PS = Pm = 1-46 x 10~^, for p2 = 0, 

P2 = 0.96x10-2, for = p5 = pAf = 0, (35) 
= 0.58 X 10~2, for pp=ps= = Pa- 



Error correction in V near D. In the presence of defects there are two modifications to error 
processes in V . First, the length scale for the minimum extension of a non-trivial error cycle shrinks. 
Second, there is a surface effect; the effective error rate for qubits in V next to the surface of defects 
is enhanced by a constant factor. 
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Figure 9: Trade-off curve for the threshold value in the presence of local and two-local errors. 
Horizontal axis: two-qubit gate error p2, vertical axis: local error rate pp + ps + Pm- 



1. Length scale for non-trivial errors: For comparison, consider a cluster cube of finite size 
2L X 2L X 2L. A non-trivial error cycle must stretch across the entire cube and thus has a weight 
of at least L. The lowest weight errors which are misinterpreted by the error correction procedure 
occur with a probability q^^'^. The total error probability incurred by such errors may therefore be 
expected to decrease exponentially fast in L, which is confirmed in numerical simulations [IT. 

In the presence of defects, the dominant sources for logical error are 1-chains that either wind 
around a defect or that begin and end in a defect and intersect a correlation surface (2-chain) in 
between; see Fig. IIUI The relevant length scales are thus the thickness (circumference) and the 
diameter of the defects. They are much smaller than L. 

Specifically, consider a defect with circumference u and length / which bounds a correlation 
surface C2, such that |{9c2}| = I. An error cycle winding around the defect has a weight of at least 
n + 8 ~ "u, and there are I such minimum weight cycles. Therefore, the probability pe{u,1) for 
affecting K{c2) by an error is, to lowest contributing order, 



In the range of validity for the above expansion in powers of q, the error is still exponentially 
suppressed in the relevant length scale u. 

2. Surface effects: As discussed above, the error level is enhanced for qubits in V near the 
surface of a defect. If the defect is primal (dual), the enhancement occurs on dual (primal) qubits. 
This effect will — if anything — lower the threshold. But there is another effect: the presence of the 
defect changes the boundary conditions. In case of a primal defect, the boundary conditions on 
the defect surface become rough for the primal lattice and smooth for the dual lattice. Dual error 
chains cannot end in a primal defect, as we noted earlier. For the dual lattice, there is excess 
syndrome available at the defect surface. This effect will — if anything — increase the threshold. Our 
intuition is that neither effect has an impact on the threshold value. The threshold should, if the 
perturbations at the boundary are not too strong, still be set by the bulk. 

We have performed numerical simulations for lattices of size L x L x 2L, where half of the 
lattice belongs to V and the other half to the defect region D. The error rate is doubled near 




(36) 
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Figure 10: Sources for logical error in the presence of defects. 
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Figure 11: Topological error correction in the presence of defects, a) Short error cycles intersect 
with the correlation surfaces and thereby jeopardize the topological error correction, b) Correlated 
errors on 5-qubits are suppressed exponentially in the qubit separation /. 



the mutual boundary of the regions and there is no remaining error in D. Simulations are feasible 
with reasonable effort up to L ~ 20. We find that finite-size effects (due to the smooth boundary 
conditions) are still noticeable up to these lattice sizes, but the intersection point of fidelity curves 
for nearby lattice sizes is slowly converging to a threshold value around that of the bulk 2.9%). 

Error correction in V near a junction between D and S. Near an S'-qubit there exist relative 
error cycles of small length, see Fig. Illb . and the topological error correction breaks down. As a 
result, the effective error on an 5-qubit is enhanced by its surrounding. To compute the effective 
error probabilities, we replace every low-weight error-chain -£(7) that results in a logical error after 
error correction by an equivalent error Esij) acting on the S-qubits. The error correction converts 
£'(7) into £'(01(7)) with 01(7) a relative 1-cycle. 'Equivalent' means that £'(01(7)) and £5(7) act 
in the same way on the stabilizer generators {-ftr^^l s G 5"} of the induced state |^'aigo)s) i-e., 
[E{ci{'~f))Es{'y), K-^ ^] = for all s £ S. The relevant correlations to check are K{c2) and K{c2) 
displayed in Fig. [TTfe .. 

It is important to note that the effective error on the S'-qubits is local. This arises because only 
error chains causing a 1-qubit error may have small length. Error chains causing a correlated error 
on the S-qubits are suppressed exponentially in the qubit separation. See Fig. Illb . 

We compute the effective error channel on the S'-qubit to first order in the error probabilities 
only. For error Model 1 the error enhancement only affects sub- leading orders of pi , 



Ti^) = (1 - p,)%] + ^ {[Xs] + [Ys] + [Zs]) . 



(37) 
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Figure 12: Shape of the defect near an S-qubit (red), a) Three-dimensional view, b) Cross-section 
through the cluster. Edges which belong to the defect are marked as "x". Black dashed lines 
connect neighbors on the cluster. Blue underlay: faces in the defect. Even/odd qubits: •/o. 



For Error Model 2 the effective error channel on an 5-qubit is not universal but depends on the 
precise shape of the defect double-tip near the 5-qubit. For our calculation we use the defect shape 
displayed in Fig.^1 The defect is one-dimensional nearest to the 5-qubit, and farther away becomes 
three-dimensional. The effective error channel then is 



= (i-Ipp-Ips- Ipm - j + y2p2 + ^ps + ^PAlj [Xs] + 
2 1 1 A 1 1 58 \ , 

■^P2 + -PS + -^PM I \Ys\ + ( -PP + -PS + ^PM + —P2 I [Zs\- 



(38) 



The individual contributions to (|38() are listed in Appendix IbI 

Error correction in S. The S-qubits are protected by the concatenated Reed-Muller code. This 
code corrects for the errors (|37|) / (|38() that remain after error correction in V and D. 

The 5-qubits are all measured in the eigenbases of ■ Then, an X- or a Y-error is equivalent 
to a Z-eiTOT with half the probability. This is easily verified for the case where X- and Y-errors 
occur with the same probability. W.l.o.g. assume the measurement basis is . Then [X] + [Y] = 

^^^^ + ^^^^ — ■^^^^(W + I-^]) — [Z]- -^^^ statement is, to leading and next-to-leading 
order in the error probability, also true when X- and ^-errors do not occur with equal probability; 
see Appendix O We may thus convert the X- and y-errors in (|37|) and (|38|) into Z^-errors. The 
corresponding error probability pz is 



2 

Pz = -pi, for error Model 1, 

76 2 4 4 
pz = —P2 + -pp + -PM + -ps, for error Model 2. 
15 3 3 3 



(39) 



The fault-tolerance threshold of the concatenated [15, 1, 3]-Reed-Muller code for independent Z- 
errors with probability pz is 1.09 x 10~^. As we discuss all Reed-Muller error correction to leading 
order only, we take the leading order estimate 

Pz,c = ~ 0.95 X 10-2. (40) 
105 
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For Error Model 1, from (|39|) and (|4flj) we obtain the threshold 



= ^ « 1.4 X 10-2, (in S). (41) 



For error Model 2 we obtain 



1 _2 

PP,c = Ps,c = PM,c = ^ ~ 0.29 X 10 , for p2 = 0, 

P2,c = ~ 0.19 X 10-2, for pp = ps=pM = 0, (in S). (42) 

532 

Pi,c = ~ 0-11 X 10"^' for pp = PS = PM = P2, 



The topological error correction in V and the Reed-Muller error correction in S run separately, 
and all that remains is to check which mechanism breaks down first. By comparison of Equation 
(j32j) with (|41|) and Equation l|35p with (|42|) . we find for both error models that the Reed-Muller 
error correction collapses first. It therefore sets the overall threshold. In Error Model 1, the critical 
error probability for local depolarizing error is 1.4 x IQ-^. In Error Model 2, for the case where 
preparation, gate, storage and measurement errors each have equal strength, the error threshold for 
the individual processes is 1.1 x 10"^. 



8 Overhead 

Denote by the number of qubits in the algorithm-specific resource state |^'aigo)Q- is also 
the number of non-Clifford one-qubit rotations in a quantum circuit realizing the algorithm "algo" . 
What is the number Nft of cluster qubits required to perform the same computation fault-tolerantly? 
The qubit overhead factor for Reed-Muller error correction is 



logiV 
logpc/p 



(43) 



where 7 = log2 15 ~ 3.91, and p is the actual and pc the critical Reed-Muller error rate. The set S 
consists of |S| = A^'Orm qubits. 

We need to determine the additional overhead due to topological error correction. We choose a 
separation r between strands of a defect loop, and a strand diameter of r/2. 



The dimensions of the cluster thus are width, depth < 3/2rNO^Mj height = 4r. The most likely 
errors then occur with a probability of exp(— K(p)r), H36() . and there are less than 3r{N0^M)'^ 
locations for them (primal+dual). We require that 3r(A^ORM)2 exp(— K(p)r) w 1, such that in the 
large N limit 

r In N. (44) 

K[p) 

Double-logarithmic corrections are omitted. The total number of qubits is given by A'^ft = width x 
depth X height < 9r^{N0^^M)'^ , such that 



iVft~iV' (logiV)3+27. 



(45) 



The overhead is polynomial. The values of the exponents in the above expression may be reduced 
in more resourceful adaptations of the presented scheme. 
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9 Discussion 



We have described a scheme for fault-tolerant cluster state universal quantum computation which 
employs topological error correction. This is possible because of a link between cluster states and 
surface codes. In addition to the topological method, we make use of a Reed-Muller quantum code 
which ensures that non-Clifford operations can be performed fault-tolerantly by local measurements. 

The error threshold is 1.4% for an ad-hoc error model with local depolarizing error and 0.11% 
for a more detailed error model with preparation-, gate-, storage- and measurement errors. We 
have not tried to optimize for either threshold value or overhead here; the foremost purpose of this 
paper is to explain the techniques. With regard to a high threshold, the obvious bottleneck is the 
Reed-Muller code. The error threshold imposed by this code on the cluster region S is — depending 
on the error model — a factor of 3 to 5 times worse than the threshold obtained from the topological 
error correction in V. To increase the threshold one may replace this code by another CSS code 
with property ^ that has a higher error threshold, provided such a code exists. Alternatively, one 
may probe the Reed-Muller code in error detection, as in magic state distillation |17j . The error 
detection threshold is 14%, which indicates that there is some room for improvement. 

Part of the investigations in this paper are numerical simulations, and we would like to comment 
on their impact on the threshold value. Numerics are encapsulated only in the threshold estimate 
for topological error correction which is much higher than the overall threshold. Our final threshold 
estimate stems from the Reed-Muller code and is analytical. 
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The appendices are relevant for Error Model 2 only. 

A Connecting sub-clusters 

The cluster C consists of a set of sub-clusters Ck, C = Ufc^fc which are prepared in sequence. Two 
successive sub-clusters Ck, C^+i have an overlap, H C^+i = C V U D. is a set of locations 
for hand-over qubits. Each graph edge (corresponding to a A(Z')-gate in cluster state creation) can 
be unambiguously assigned to one sub-cluster. The set of edges ending in one vertex either belongs 
to one or to two sub-clusters. In the latter case, the vertex is the location for a hand-over qubit. 

If k is odd, the cluster state creation procedure on is the sequence 1 — > 2 — > 3 ^ 4. If 
k is even, the sequence is 3 — > 4 — > 1 — > 2. As shown in Fig. 113b . connecting the sub-clusters 
proceeds smoothly. Consider, for example, the sub-cluster Ci. With the exception of a subset of 
the hand-over qubits to C2, the qubits in Ci are prepared at t = 0, entangled in steps 1 to 4, 
and measured at t = 5, 6. Specifically, the V- and £)-qubits are measured at time t = 5, while 
the S'-qubits are measured at time t = 6 with measurement bases adapted according to previous 
measurement outcomes. The sub-clusters Ck are chosen such that there is no temporal order among 
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Figure 13: Connecting sub-clusters, a) Elementary cell of L. The old sub-cluster is displayed in 
blue, the new in black; "o"; hand-over qubits, qubits within one sub-cluster, b) Temporal order 
of operations. "P": |+)-preparation, "E1-E4": steps of parallel A(Z)-gates, Ml, M2: measurement. 



the measurements on qubits within 5" H C^, for all k. Then, the 5'-qubits wait for one time step in 
which a storage error may occur. 

The hand-over qubits stay in the computation no longer than the other V- and D-qubits. They 
form two subsets, }i\ = A\J B\ see Fig. I13t i. The qubits in A are prepared at t = 0, entangled in 
time steps 1 to 4 and measured in at t = 5. They cause no change in the error model. 

The qubits in B are not acted upon by a gate until time step 2 (since the potential interaction 
partner of step 1 isn't there yet), so they are prepared at time t = 1. The final interaction involving 
the i?-qubits is in step 5, and they are measured in step 6. Between preparation and measurement, 
the -B-qubits are in the computation for four time steps in each of which they are acted upon by 
a gate. No additional storage error occurs. There is one modification due to the i?-qubits. The 
temporal order of A(Z)-gates involving the qubits 6 G is changed. As a result, the correlated 
errors on the edge qubits of the faces {02} = h are not among pairs of opposite edge qubits but 
among pairs of neighboring edge qubits. So, the error rates q\ and ^2 in are unchanged, but 92 
characterizes a slightly different process. 

We expect this to be a minor effect. The overall threshold is still set by the threshold for Reed- 
Muller error-correction, which is some five times smaller than the simulated threshold for topological 
error correction. 



B Effective error channel on the /S-qubits 

The effective error on an S'-qubit stems from the S'-qubit itself and its immediate surrounding shown 
in Fig. and from the two edge-qubits in the one-dimensional section of the defect, which are not 
protected by any syndrome. Of the latter each contributes an error 

/4 2 2 \ 

£^defect = ( 5^*2 + -^VS + l^VM \ \Xs\. (46) 

The preparation error does not contribute, because the corresponding Z-error on the defect qubit 
is absorbed in the Z-measurement. 

The effective error of the center qubit s £ S stems from operations that act on s directly, from 
X- or y-errors propagated to s by the A(Z)-gates and from short nontrivial cycles. Specifically, 
there are four non-trivial cycles of length 3. One of them is denoted as E{ci) in Fig. 1141 Because of 
the correlations in the forward-propagated errors these cycles have weight 2 and cause inconclusive 
syndrome at lowest order in the error probability. There are further error cycles of length 3, such as 
E{ci) in Fig 1141 But they have weight 3 even for Error Model 2 and do not contribute to the lowest 
order error channel. We perform a count including all error sources in the cluster region displayed 
in Fig. [Til right. There is one convention that enters into the count. Namely, the error ZnZiyZ^Zfi 
(see Fig. I14j) is a non-trivial error cycle such that the errors ZaZf, and ZcZ^, have the same weight 
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defect 




Figure 14: Lattice around an S'-qubit. Left: Location of the S'-qubit (red) and the defect (gray). 
E(ci) is an error cycle of length 3 and weight 3 for Error Models 1 and 2. Right: Detail. E{ci) 
is a non-trivial error cycle of length 3 which has weight 2 in error Model 2. Red: S'-qubit, black: 
lattice C for topological error correction. The cluster edges (blue) correspond to A(Z)-gates whose 
temporal order is as indicated. 



and the same syndrome but different effect on the computation, Xg vs. Ig- When the corresponding 
syndrome occurs, we assert that a logical X-error occurred and correct for it. We obtain 

(211\ /211 58 \ 

■^P2 + ^PS + ^PM U[Xs] + [Ys] )+i-pp + -PS + ^PM + -^P2 j[Zs]. (47) 

The sources (|l6|) and (|T7|) combined, -^central + 2 x E'defect, lead to the local error channel (|38|l . 

C Conversion of X- and y-errors on S'-qubits 

Here we show that an X- or Y-error on an individual S'-qubit s with probability p is equivalent to a Z- 
error on that qubit with probability p/2, p[Xs] = p[Ys] — p/2 [Zg]- The qubit s may be measured in 
the eigenbasis of or of . W.l.o.g. assume the qubit s is measured in the eigenbasis of . 

Then, Xg = ^"J^'' ~ ^°v^^° " '^^^ ^-error is equivalent to a coherent Z-error and we need to 

check whether the coherences matter. More generally, for the described scenario with the subsequent 
measurement a probabilistic local error channel (1 — px — Py — Pz)[Is] +Px[Xs] + pY[Yg] +pz[Zg] 
is equivalent to a channel with coherent errors 

P — > {l-q)p + qZgpZg + iq [pZg - Zgp) , (48) 

with 

. Px+PY ~ Px-PY 
q = Pz + ^ , q = ^ • (49) 

Now assume that all S-qubits are affected individually by the error channel 1)48(1 . The Reed-Muller 
error correction at successive levels maps these channels to channels of the same form, one coding 
level higher up. The parameters qi, qi at coding level I obey recursion relations which, up to fourth 
order, read 

qi+l = W5qf {I- qiY^ +35 qf {I- qiY^ + 1260 qf (1- qi)^^ +630 qf {I- qi)^\ 
qi+l = 70i{l-qi)^^ -16^0 iqi {I- qif\ 

If we compare ()5U() to the recursion relation of q for probabilistic Z-error, a deviation first shows up 
at fourth order. The discussion of Reed-Muller error correction in this paper is confined to leading 
order. The leading order result for the threshold, qc = 1/105, is not affected by the coherences 
in the error ()48|) . Probabilistic X-and Y-errors influence the threshold by contributing half their 
weight to the probability of an effective Z-error; see (|49p . 
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